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We demonstrate electromagnetically induced transparency (EIT) in a four-level cascade system where the 
upper level is a Rydberg state. The observed spectral features are sub-Doppler and can be enhanced due to 
the compensation of Doppler shifts with AC Stark shifts. A theoretical description of the system is developed 
which agrees well with the experimental results and an expression for the optimum parameters is derived. © 
2013 Optical Society of America 
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Non-linear optical effects in ensembles of Rydberg 
atoms is a topic of burgeoning interest [1]. The large 
polarizability of the Rydberg states leads to very large 
DC electro-optic effects [2] and the strong, long-range 
interactions between neighboring atoms have been ex- 
ploited to demonstrate cooperative enhancement of non- 
linearities [3,4] and single photon generation [5]. Ry- 
dberg non-linear optics offers the prospect of non-local 
effects [6] , photon-photon interactions [7] and single pho- 
ton subtraction [8]. 

Multi-photon excitation schemes are of particular in- 
terest because they can eliminate motional dephasing of 
dark state polaritons [9] and allow Rydberg states to be 
accessed using inexpensive and convenient diode laser 
systems. Recent experiments have demonstrated Ryd- 
berg spectroscopy based on population transfer [10, 11], 
and similar coherent processes have been studied theo- 
retically [12]. 

Here we study, experimentally and theoretically, co- 
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Fig. 1. a) Experimental setup, b) Bare state atomic 
energy- levels, c) Dressed state energy levels; the dressing 
laser acts as a perturbation to the system. 
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Fig. 2. Experimental weak probe transmission spectrum 
for P p = 20 nW, P d = 500 /xW and P c = 200 mW. 



herent three-photon electromagnetically induced trans- 
parency (EIT) [13] involving Rydberg states in Cesium 
vapor. A schematic of the experimental setup and the 
atomic energy level scheme are shown in Fig. 1 (a) and 
(b) respectively. States |1), |2), |3) and |4) represent the 
6Si/ 2 , 6P3/2, 7Si/2 and 26P 3 / 2 states respectively. A 
weak 852 nm probe beam coupling states |1) and |2), is 
focused through the center of a 2.4 mm Cs vapor cell and 
then detected at a photodiode. The counterpropagating 
1470 nm dressing beam is stabilized to the |2) — » |3) tran- 
sition using excited-state polarization spectroscopy [14]. 
The 790 nm coupling beam also counterpropagates with 
the probe and excites the |3) — >> |4) transition. We scan 
over either the |1) — )• |2) or |3) — >• |4) transition whilst the 
other two lasers are on resonance. The 1/e 2 radii of the 
probe, dressing and coupling beams are 14 /im, 30 /im 
and 20 /im, the wavevectors are k p , kd and /c c , the beam 
powers are P p , Pd and P c and the Rabi frequencies are 
£7 P , and £l c respectively. The temperature of the cell 
is maintained at 50° C giving a number density, Ao, of 
2.0 x 10 11 cm- 3 . 

A typical spectrum obtained when the probe fre- 
quency is scanned is shown in Fig. 2. A narrow spec- 



1 




-50 -25 25 50 -50 -25 25 50 -50 -25 25 50 

Coupling Detuning (MHz) 



Fig. 3. (color online) Top: Absorption coefficient per velocity classes from the four-level model for different values 
of Qd/Q c , (a) 0.6 (b) 1.2 and (c) 2.5. At the optimum ratio of 1.2 (calculated from Eq. (7)), the bright state is 
insensitive to the atomic velocity. Bottom: Typical experimental spectra with fits from the full theoretical model for 
different values of dressing laser power Pa, (d) 17.4 /iW, (e) 79.7 pW and (f) 345.2 pW. Experimental parameters 
fixed at P p = 20 nW, P c = 41.4 mW and A p = A d = MHz. 



tral feature is observed within the D2 Doppler absorp- 
tion lineshape with enhanced absorption on resonance 
and enhanced transmission either side of resonance. This 
lineshape can be understood by considering the dressed 
state picture shown in Fig. 1 (c) where the system now 
contains two cascade EIT systems [15]. The result is two 
dark states which overlap in frequency space and inter- 
fere constructively to give enhanced or electromagneti- 
cally induced absorption (EIA) [16]. When the coupling 
laser is scanned instead, a similar lineshape is observed 
but on a fiat transmission background as shown in Fig. 

3 (d-f). 

To model the system we consider the four-level sys- 
tem shown in Fig. 1 (b). Using semi-classical theory and 
including the Doppler effect [19], the complex suscepti- 
bility of the system per velocity class in the weak probe 
approximation is given by 
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where d 2 \ is the probe dipole matrix element and N(v)dv 
is the atomic density per velocity class v. This Gaussian 
distribution is written in terms of the most probable ve- 
locity u = ^/2k^T /m and given by 
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The density matrix elements p 2 i, P31 and p^i represent 
the coherences between the ground state |1) and excited 
states 1 2), |3) and |4) respectively. The ratio psi/ ' P21 be- 
tween the coherences is 
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where the coherence ratio p^\j pz\ is 
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A c is the coupling beam detuning from resonance and 7 C 
is the dephasing of the Rydberg state due to transit time 
broadening and laser intensity variations. T2 and T3 are 
the decay rates of states |2) and |3) respectively. The 
total complex susceptibility Xtot is given by integrating 
Eq. (1) over all velocity classes. The total absorption 
coefficient is given by 



«tot = fe p Im[xtot]. 



(5) 



Fig. 3 (a), (b) and (c) show the absorption coefficients 
per velocity class from the simple four-level theoretical 
model for Q^/Q c = 0.6, 1.2 and 2.5 respectively as a 
function of coupling field detuning A c and atomic veloc- 
ity v. We can see that, for different values of ^d/^c, the 
gradient of the line defining the frequency of maximum 
absorption rotates from positive (Fig. 3 (a)) to infinity 
(Fig. 3 (b)) to negative (Fig. 3 (c)). Significantly, at the 
optimum ratio of Qd/Q c , the resonance becomes veloc- 
ity insensitive for A c = MHz. As a result, more atoms 
contribute to the signal at the optimum ratio and the ab- 
sorptive resonance reaches its maximum magnitude, as 
shown in Fig. 3 (d-f). This occurs because the Doppler 
shifts of the atoms are compensated by AC Stark shifts 
[18]. The theoretical fit curves in Fig. 3 (d-f) are calcu- 
lated by averaging over all magnetic sublevels [17]. 

As increases, the energy separation between the 
two dressed states |a + ) and |a _ ) increases, resulting in 
the increase in frequency separation between the two EIT 
resonances as shown by Fig. 4. This splitting scales with 
the square root of power and has a minimum splitting 
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of approximately 26 MHz which is determined by the 
wavevectors and Rabi frequencies of the laser fields [19]. 

To understand the optimum ratio of Q^/Q c , we con- 
sider the Hamiltonian of the bare state system H given 
by 
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where Ai ph = -k p v, A 2ph = ~(k p - k d )v and A 3ph = 
— (fc p — fed — k c )v. To calculate the eigenenergy of the 
dressed state, we diagonalize the Hamiltonian contain- 
ing the three upper states and consider the eigenenergy 
around v = m/s. We can then expand the character- 
istic equation in powers of v and neglect higher order 
terms. The resulting condition for a velocity insensitive 
bright state is 
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Eq. (7) sets the condition in which the absorption 
reaches its maximum level as many velocity classes now 
contribute to the total absorption at A c = MHz. Note 
that this condition predicts the optimum condition for a 
single four-level system. However in the real system and 
our full model, we have many different four-level sys- 
tems distributed over the magnetic sublevels. The mag- 
nitude of the spectral feature as a function of P d for P c 
= 41.4 mW is shown in Fig 4. The spectral feature has a 
maximum magnitude at P d = 79.7 fiW. At this point, for 
the strongest transition, the ratio ftd/ft c = 1-1 which is 
close to the value of ftd/ft c = 1-2 predicted by the simple 
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Fig. 4. (color online) Maximum change in transmission 
(red squares) and transparency splitting (blue circles) as 
a function of The solid (red) line is the theoretical 
calculation from the full model. The dashed (blue) line 
is a scaling fit with power. The optimum ratio of ftd/ft c 
occurs at P d = 79.7 /iW. 



four- level analysis in Eq. (7). The difference between the 
actual value and the value predicted by the simple model 
is due to the inclusion of transitions from all magnetic 
sublevels. 

In summary, we have demonstrated Doppler- 
compensated Rydberg EIT. We show that there is an 
optimum ratio of Rabi frequencies for the upper two 
steps which leads to an enhanced sub-Doppler feature. 
This work could provide an important step towards 
single-photon non-linearities in the telecoms wavelength 
range. 
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